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Abstract We analyzed 186 binary pulsars (BPSRs) in 
the magnetic field versus spin period (B-P) diagram, 
where their relations to the millisecond pulsars (MSPs) 
can be clearly shown. Generally, both BPSRs and 
MSPs are believed to be recycled and spun-up in binary 
accreting phases, and evolved below the spin-up line 
setting by the Eddington accretion rate (M ~10 18 <?/s). 
It is noticed that most BPSRs are distributed around 
the spin-up line with mass accretion rate M = 10 le g/s 
and almost all MSP samples lie above the spin-up line 
with M ~ 10 15 <7/s. Thus, we calculate that a minimum 
accretion rate (M ~ 10 15 g/s) is required for the MSP 
formation, and physical reasons for this are proposed. 
In the B-P diagram, the positions of BPSRs and their 
relations to the binary parameters, such as the com- 
panion mass, orbital period and eccentricity, are illus- 
trated and discussed. In addition, for the seven BPSRs 
located above the limit spin-up line, possible causes are 
suggested. 



Keywords pulsars: 
neutron 



general; binaries: close; stars: 



1 Introduction 

The pulsars (PSRs) in the binary systems are al- 
ways convenient for us to determine their parameters 
such as the mass, origin and evolution. The major- 
ity of the binary pulsars (BPSRs) are millisecond pul- 
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sars (MS Ps) with the spin period less than 20 mil- 
liseconds ( Stairs! 12004 ). Based on the data from the 
ATNF Pulsar Catalogue in November 2012, it has 
been found 2008 PSRs including 221 MSPs (136 bi- 
nary PSRs versus 85 isolated ones). The compan- 
ions of 186 binary PSRs are white dwarfs (170), neu- 
tron stars (9), main sequence stars (4) and planets 
(3), respectively. Much progress has been achieved 
in understanding the form ation and evolution of the 
binary pulsars (BPSRs) (IBhattacharva et al.l Il995 



Stairj2004tJManchesteiJl2004[|Lorimer fc Duncanll2008 



TaurislboialTauris. Langer fe Kramer II2012J ) 

The distributions of the magnetic field (B) and 
spin period (P) are both almost bimodal with a di- 
chotomy between the normal pulsars (PSRs) (B ~ 
W n - la G, P ~ 0.1s - 10s) formed directly by the 
supernova explosions and MSPs (B ~ 10 8_9 G, P ~ 
1.4ms — 20ws) that are recycled in the binary sys- 
tems ( Aloar et al. 1982t Radhakrishnan fc Srinivasan 



1982 [lLorimer fc Duncar]|2008l : lKaspill2010l : IWang et al 



20111 ). In the binary systems, with the accreting mat- 
ter of ~ 0.1 — 0.2M Q from their companions, the neu- 
tron stars (NSs) can be spun-up to several millisec- 
onds, while their magnetic fields decay to ~ 10 8 ~ 9 G 



( Bhattacharva fc van den Heuvelll99ll:JBhattacharva et al.l 



1995 



2005 



van den Heuvel fc Bitzarakill 995: M elatos fc Pavnel 
Zhang fc Koiimall200rf: IWang et al.ll201ll) . With 



the 3D simulations of magnctohydrodynamics (MHD), 
the magnetic fi eld affected by the accretion fl ow is de- 
picted clearly ( Kulkarni fc Romanova! 120 08) . When 
the accretion phase is finished, the radio emissions 
of the fast rotating NSs can be detected as MSPs 
( Alpar et al. 1982t Radhakrishnan fc Srinivasan! Il982 : 



Tauris. Langer fc Kramer ll2012l ). 

Evaporation of the companion star by the pulsar 
radiation explain explains why som e MSPs are single 
( Kluzniak et alJl98alHessels ll2008h . This is more com- 
plicated in globular clusters: a companion can be dis- 



rupted by encounter with another star. Direct evidence 
supporting the MSP spin-up model has been found, for 
instance, the accreting X-ray MSP SAX J1808. 4-3658 
(B ~ 10 8 G, P = 2A 9ms) discovered in the low mas s 
X-ray binary system (jWiinands fc van den Kilslll998| ). 
and the double pulsar PSR 0737-3039A/ B (a MSP and 
a normal PSR in the binary system) (IBurgav et al 



20031 Ivan den HeuvellliooS ILvne et aill2004t) 

BPSRs are the most accurate clocks in the unverse 
that can be ex ploited to form a MSP detect gravi- 
tatio n al waves (iManchesterl l2004t lLorimer fe Duncan 
2008J iHobbs et all |2010|) . BPSRs are the most ac- 
curate clocks in the unverse In addition, the BP- 
SRs are significant systems for testing general relativ- 
ity effects, e.g., the gravitational-radiation-induced or- 
bit shrinking, gravitatio nal red-shift and Shapiro de- 
lay ( Kramer et al.ll2006l) . based on which the NS mass 
can be measured in double neutron star systems us- 
ing the periastron advance and with the high precision 
( Zhang et ai1l201ll) 

This paper mainly focuses on the analysis of 186 BP- 
SRs in the B-P diagram with the data from ATNF pul- 
sar catalogue (Manchester et al. 2005). To understand 
the properties of these pulsars, we pay attention to their 
B-P positions relative to the spin-up lines, in relation 
to the companion mass, eccentricity and orbit period. 
The structure of the paper is organized as follows: 
Sec. 2 presents the detail study of the spin up lines un- 
der the different accretion rates. Here we also give 
the implications of the minimum accretion rate on the 
millisecond pulsar formation. In Sec. 3, the distribu- 
tion characteristics of BPSPs in the B-P diagram as a 
function of the binary parameters and the evolution- 
ary scenario of binary pulsar are discussed. In Sec. 4, 
we suggest the causes why seven BPSRs lie above the 
Eddington limited spin-up line. The discussions and 
conclusions are given in Sec. 5. 



2 Pulsars in B-P diagram and the evolutionary 
implications 

2.1 The Spin- up line 

If a NS, which has a symmetry axis of the magnetic 
field aligned with the rotation axis, spins up to the 
Kepler orbital period at the magnetospherc radius, 
a spin-up line will be obtained. It is also called 
the equilibrium spin period line expressed the rela- 
tions hip between the magnetic field and the spin pe- 
riod dGhosh fc Lamblll977J IShamro fc Teukolskvlll983l: 



Bhattacharva fc van den Heuvellll99ll: iGhosh fe Lamb 
1992tlBurderi. D'Antona fe Burgavl2002l:lFrank. King fc 



2002t ICamenzind M.1 120071 ). However, the magneto- 
sphere radius depends on many factors, such as the 
structure of the accretion disk and mass transfer form, 



etc (jTauris. Langer fc Kramer Il201ll 120121 ). 

The equation of the spin-up line ca n be written as 
( Bhattacharva fc van den Heuvellll991 ): 
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where P eq is the "equilibrium" spin period, B 9 is the 
dipole magnetic field in units of 10 9 G, M is the accre- 
tion rate, M^dd is the " Eddington- limit" accretion rate 
(10 18 g/s), i? 6 is the stellar radius in units of 10 6 cm and 
M is the NS mass. 

The spin period (P) and the spin-down rate (P) 
of a PSR can be combined to yield an estimate of the 
strength of the dipole magn etic field ([Manchester fc Tavlor | 
1977t ILvne fc Smith! 120061) : 
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where / is the moment of inertial and R is the radius 
of the PSR. The moment of inertia is / = (2/5)MP 2 . 
To evaluate the influence on the divided magnetic field 
by the choice of NS mass and radius, we rewrite Eq.Q 
in terms of M and R: 



B = 3.4 x 10 ia ( 
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Eqs. (fiT3"l) show us that the NS mass (M), radius (R) 
and accretion rate (M) influence on the position of the 
spin-up line in the B-P diagram, and both NS masses 
and radius have influences on the derived strength of 
the dipole magnetic field. Mostly one takes a standard 
pulsar parameter (M = IAM@, R = I0 6 cm) in the 
three equations. However, some measured NS masses 
are larger or smaller than 1.4Mq. In binary systems, a 
mean value of M = 1.46±0.30 M^ has be en derived for 



the 61 measured NS masses (jZhang et al.l 120111 ). The 
NS radius has not yet been measured as accurate as its 
mass, and is usually esti mated in the range 10 ~ 20km 



by means of the models (Miller. Lamb fc Psaltisl 



Miller 



2007). 



1998 



2002t lLarrimer fc Prakasbl 12004 lliessels etal 



2006 : lHaensel. Potelin fc Yakovlevl 120071: I Zhang et al 



The spin-up line equation Eq.(l) and magnetic field 
equation Eq.(3) show the relationships between the 
PSR magnetic field and NS radius as B ex R~ 3 and 



R|me 



Li? , respectively. If a NS radius changes from 10 

I 

To 20 km, then the spin-up line will shift down by 



Log(8)=0.9, whereas the PSR position will shift down 
by Log(4)=0.6, in Log(B) axis of the B-P diagram. 
That is to say, the spin- up line will shift down ~ 0.3 
relative to the PSR position in the Log(B) axis of the 
B-P diagram. 

Compared with the affection of the radius on the 
spin-up line and PSR position in B-P diagram, the vari- 
ation of accretion rate will give much more influence on 
the spin-up line. From the observations, the accretion 
rates of the X-ray NSs span th ree orders of magnitudes, 
e.g., f rom 10 15 <7/s to 10 1 8 a/s (IHasinger fc van der Klis 
1989t Ivan der Klisll200ol lLamb fe Yull2004l ). and these 
variations will give rise to a considerately changes in 
the positions of the spin-up line in the B-P diagram, 
which are shown in FigO]for accretion rates: 10 15 g/s, 
10 16 g/s, 10 1 V s and W ls g/s. 

The distribution of the isolated and binary PSRs are 
also shown in Fig[TJ The majority of the MSPs are dis- 
tributed around the spin-up line with M = 10 16 g/s. It 
seems likely that most MSPs may be born at the posi- 
tions near the spin-up line with M = 10 17 g/s, which is 
the generally observed accretion r ate for the X-ray NSs 
in the Low-Mass X-ray Binaries (Bradt fc McClintock 



minimum mass accretion rate for the MSP formation 
can be estimated as: 



19831: IHasinger fc van der Klislll989Hyan der Klisl l2000l: 



Liu, van Paradiis fc van den Heuvell l2007h . This fact 



was also pointed out a decade ago by Lundgren et al. 
(1996) and Hansen and Phinncy (1998). 

2.2 Implications of the minimum accretion rate for 
MSP formation 

It is noticed that only three MSPs, J0514-4002A, J1801- 
3210 and J1518+4904, lie below the spin-up line with 
M = 10 15 g/s. With the data from ATNF and Eq.flU), 
their positions correspond to the mean accretion rates 
0.6 x 10 15 g/s,0.8 x 10 15 g/ s and 0.9 x 10 15 g/s, re- 
spectively. Possibly, their spin periods evolved to the 
present values due to the spin-down by of the electro- 
magnetic emission. This suggest that there exists a 
critical minimum accretion rate M = 10 15 g/s for a 
MSP formation, and the arguments for this view are 
described in the following. 

At first, the real ages of MSPs have not yet been 
satisfactorily determined, although their charactcris- 
tic ages can be a little longer tha n the Hubble age 
( Camilo. Thorsett fc Kulkarnilfl994l) . During the MSP 
spin-up evolution, the minimum amount of accreted 
mass required for a MSP to reach the spin period of 
several milliseconds is about AM rr ' ■> 0.1 - 0.2M Q 



( Wang et all 1201 it iTauris. Langer fc Kramer I |2012J), 



and the time of spin-up should be less than the Hubble 
time of about tu = 10 10 yrs of a MSP. Thus, a critical 



M r 



AM cr /t H ~ Wg/s 



(4) 



If the accretion rate is much lower than this critical 
value, the NS will have no chance to accrete sufficient 
mass to be spun-up to a period of a few milliseconds 
during the evolution. Based on the above argument, a 
critical minimum accretion rate ~ 10 15 g/s is suggested 
as a necessary condition for the MSP formation. 



3 Binary Pulsars in the B-P Diagram 

The evolution of the binary parameters of BPSRs, such 
as the companion mass, eccentricity and orbital period, 
are topics of interest. However, the BPSRs in GC can 
not help us to understand their properties ver y well 
( Hessels Il2008t ITauris. Langer fc Kramer 112012 ) since 
the stellar density of the cluster is much higher than 
that of the Galactic disk which increase the collision 
and capture of the stars. In this part, we ignore the 75 
BPSRs in the GC, and study the distribution of 103 
BPSRs in the Galaxy disk (GD) that have spin period 
P and magnetic field B in the B-P diagram. 



3.1 BPSRs with the Different Parameter Conditions 

In Fig. [5^,, we notice that the BPSRs with the small 
and inter-mediate companion mass, M c < 0.1 Mq and 
0.1M Q < M c < 0.45M© as defined by Stairs (2004) 
occupy about 12.6% and 61.2% in the total BPSRs, 
respectively. In FigJ^b, the percentages of the BPSRs 
with the low eccentricity (e < 0.1), mediate eccentricity 
(0.1 < e < 0.5) and high eccentricity (e > 0.5) occupy 
85.1%, 6.9% and 8.0% of the total, respectively, which 
indicates that the eccentricity of BPSR has deceased in 
the evolution process. In FigJ^fc, most of the BPSRs 
with the orbital period P or b < lOOd share the spin pe- 
riod less than ~ 100ms. Then those with long orbital 
period P or \, > lOOd are distributed from ~ 30ms to 
~ls. 

3.2 Evolution Scenario of Pulsar in Binary System 
Ther e are many cases for the evolutions of binary sys- 



tems ( Lipunov fc Postnovll984l : lTauris. Langer fc Kramer | 



20111 120121 ) . The companion mass p lays an important 

role in the evolution of binar y system ( Taam fc van den Heuvel^ 
1986tlShibazaki et al.lll989l ). The star with a mass over 
about 8 solar masses experiences a supernova explo- 
sion firstly, leaving a NS with a high magnetic field of 
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Fig. 2 B-P diagrams plotted with various parameter conditions of 103 BPSRs, as shown in the figures, e.g. companion 
mass (a), eccentricity (b) and orbital period (c). Ab out 86.4% of the BPSRs have low m agnetic fields B < 10 10 ' 5 G, which 
is considered as the threshold for the normal PSRs (|Hartman, Portegies & Verbuntlll997r ). 
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Fig. 1 The magnetic-field versus spin-period diagram for 
2008 PSRs (data from ATNF pulsar catalogue as at Novem- 
ber 2012). 221 MSPs (P < 20ms) are shown, including the 
isolated ones (85) and those in the binary systems (136). 
The solid, dash, short dash and dash dot-dot lines stand for 
the spin-up lines with the different accretion rates 10 18 <;/s, 
10 17 g/s, 10 16 gr/s and 10 15 g/s, respectively, as shown in Eq. 
(1). The dash dot line is the death line by the theoreti- 
cal model of deceasing the radio emission of pulsar (Bhat- 
tacharya & van den Heuvel 1991; Ruderman & Sutherland 
1975). The isolated and BPSRs are labeled by dots and 
stars, respectively. 



~ 10 12 13 G and spin period of > 10 milliseconds, like 
the Crab pulsar {B ~ 10 12 G, P - 30ms). Such a new 
born pulsar in the B-P diagram lies in the region A0 
in FiglSl from where the pulsar starts to evolve to Al 
on account of the electromagnetic emission with a more 
or less constant magnetic field within the time about 

- ioV- 

The life time of NS spin-up will be inversely re- 
lated to the mass of its companion, since heavier 
stars have shorter life spans, which can be an ap- 
proximated by the formula of the main sequence age 
T - lO lo (2/r)(A//M )" 2 ' 5 . With the accretion mass 
provided by the companion, the NS will be spun up. 
The e volution can be classif i ed in the following three 
types (jShibazaki et al.l Il989t iTaam &: van den Heuvel 
19861) . 



(1) High Mass X-ray Binary (HMXB) case: the com- 
panion mass is about 10-25M Q . The binary system 
will be the PSR+massive (PSR+NS) star if the sys- 
tem gets the accretion mass ~ 0.0001M Q — 0.001M Q 
(~ 0.01M Q ), while the evolution path might be A1-A2 
(A1-A2-A3). The accretion evolution will stop at the 
spin-up line, where the stellar spin frequency equals the 
orbital Kcplerian frequency at the edge of the magne- 
tosphere. 

(2) Intermediate Mass X-ray Binary (IMXB) case: if 
the companion mass is in the range of 3-1OM0, the NS 
spin-up evolution will continue the procedure of HMXB 
case following up the route A1-A2-A3 or even to posi- 
tion A4 along the spin-up line. The NS will accrete 
about - 0.001M Q -0.1M Q to form a NS+ WD (heavier) 
systems, with the magnetic field B ~ 10 9-10 G and the 
spin period P ~ 20 — 100ms. 

(3) Low Mass X-ray Binary (LMXB) case: if the 
companion mass ranges at ~ 0.8 — 3M Q , the NS will 
evolve to the position A4 with the magnetic field B ~ 
10 8_9 G and the spin period P < 20ms, while the NS 
accretes > - 0.1 - 0.2M W to form a MSP+WD(lighter) 



binary system (jZhang et al.ll2011l) 



19861) . 



From the above mentioned scenarios of BPSR cvo- we think that it has little possibility to accrete the mat- 
lutions, the companion mass and accretion rate might ter from its companion star, therefore it is inferred that 
lay important roles in the evolution of binary system th is pulsar is a non-recycled BPSR, born at the position 
Wiierslll997[IShibazaki et al.lll989tlTaam fc van den Heuvedfr ar AO. 

The common characteristics of these four BPSRs are: 
they all have high eccentricities, long orbital periods, 
massive companions in the main sequence stage. In 
all cases the pulsars are new born ones with spin-down 
ones of less than a million years,; they are not (yet) 
accreting, and are in the "non-recycled" phase. These 
four BPSRs are on the way of evolution from AO to 
Al as shown in FigJSl where the first born PSRs are 
experiencing the spin-down with little accretion mass. 



4 Binary pulsars above he spin-up line 

In the binary systems, generally, the PSRs will evolve 
below the spin-up line after accreting a sufficient 
amount of mass from their companions. However, 
seven BPSRs lie above the Eddington spin-up line 
(M = W 18 g/s) with the magnetic fields > 10 115 G, as 
shown in Figj4] Their parameters are listed in Tab. 1. 
According to their companion masses, these seven BP- 
SRs are divided them into two types: the ones with the 
massive stars (M > 4.0M Q ) and the others with the 
degenerate stars (NS or WD). 

4.1 Type-I: BPSRs with Massive Companions 

Type-I includes four binary pulsars. Three of them 
(No. 1,3 a nd 4 in Figj i p are with B or Be sta r com- 
panions (Kaspi et al. Il994t Manchester et al.1 2001 : 
Stairs. Manchester fc Lvne! l200lt Istairs et all 12003 : 



Wang et al.ll2004h . PSR J1638-4725 (No.2 in FiggJ has 
a massive star as its companion (M c > 8.08A/ Q ) which 
is in the main sequence stage. Its radio emission dis- 
appears near the periastron because of the absorption 
and scattering by the dense stellar environment. Dur- 
ing the radio-quiet phase, X-rays that are produced by 
the accretion matter o nto the NS magne tosphere are 
likely to be detectable (JMcLaughlinl 120041 ) . By consid- 
ering the parameters of this pulsar, such as the strong 
magnetic field (B ~ 1.93 x 10 12 G), not long charac- 
teristic age (r = 2.53Mj/r), very long orbital period 
(Porb = 1940.9c?) and very high eccentricity (e = 0.955), 
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4.2 Type-II: BPSRs with Degenerate Star 
Companions 

J1906+0746 (No.5 in Figg} is a non-recycled PSR 
(1.25M Q ) in a double neutron star system. Its com- 
panion is a recycled PSR with the mass of 1.37M Q 
( Lorimer fc Duncanll2008l ) . The short characteristic age 
of J1906+0746 (r = 113000yr) indicates that it is a re- 
cently formed young PSR after the core collapse, which 
should be the reason why its B-P position lies above 
the spin- up line (10 18 g/s). 

B1820-11 (No. 6 in Fi glU) possesses a slightly massive 
WD (M c = 0.78 Mfn) (jPortegies. Simon fc Yungelson 
1999t IStaira I2004T ) as its companion. From its param- 
eters (t = 3.22Myr, B = 6.29 x 10 n G and P = 
279.829ms), we suggest that it is a young recycled PSR. 
The evolution history of this PSR can be understood in 
this way: the initial magnetic field of NS can be as high 
as B ~ 10 13 G, and it evolves from the position Al to A2 
after accreting the companion mass about ~ 0.001M©, 
while one to two magnitude orders of magnetic field 
has been deducted. Due to its high eccentricity and 
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3 Evolution tracks of binary pulsars in B-P diagram. 
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Fig. 4 B-P diagram for 103 BPSRs with the spin period 
and magnetic field in GD. The plots share the similar mean- 
ings to those of Fig.l. 



large orbital period, B1820-11 might not get enough 
accreting mass from its companion. In Eq.([T]) it is as- 
sumed that the spinning axis is parallel to the magnetic 
axis, whereas in Eq.@ the perpendicular condition is 
assumed. Therefore the position of the spin-up line will 
be fuzzy because the angle of the spinning axis to the 
magnetic axis is unknown. The position of PSR B1820- 
11 might be below the spin- up line. 

J1141-6545 (No.7 in Figg} is a PSR of mass 1.3 
M Q with an optical W D (1.02M^) as its com panion 
in the binary system ( Antoniadis et al.l 12011 ). With 
the short orbital period {P or bit = 0.1977(f) and slightly 
low eccentricity (e = 0.1719), it can be derived that 
this PSR acquired the accretion mass easily and fol- 
lowed up the recycled process. Its position upon the 
spin-up line can be explained by two possible rea- 
sons: the first one is the same reason as for B1820- 
11, and the second one is that the PSR may have 
experienced the acc reting induced collapse (AIC) of a 
white dwarf (WD) (Ivan den Heuvel &: Bitzarakilll995 ; 
van den Heuvel 200 



It is noted that the spin- up line with Eddington rate 
in B-P diagram is also influenced by the spin-up torque 
and radiation pressure i n disk, which has been stud- 
ied by some researchers (IKulkarni fc Romanoval 12008 ; 
Ghosh fc Lamb3ll992t [Lamb fc Yull2004l) . 



5 Discussions and Conclusions 

In this paper, we investigate the 186 BPSRs (136 MSPs) 
in the B-P diagram. The origin and evolution informa- 
tion of their progenitors can be inferred by analyzing 
their positions relative to the spin-up lines in the B-P 
diagram. The following results have been obtained: 

(1) The observations of X-ray NSs in binary systems 
show that M -- 10 17 q/s is a conventionally detected 
. der Klisll2000tlLiu. van Paradiis 



mass accretion rate 



20071) . In the B-P diagram, most BPSRs are distributed 
around the spin-up line for M = 10 16 g/s, which means 
that most BPSRs should be born at the positions near 
the spin- up line (M = 10 17 g/s), from where the PSRs 
experience the spin-down time of ~ 10 8_9 j/r after the 
onset of the radio PSR emissions. 

(2) To form a MSP, the mass of about at least 
~ 0.1-0.2 M0 should be accreted from the companion 
( Wang et al1l201ll Izhang et al.ll201ll ). Thus, within a 
Hubble time, the required minimum accretion rate is 
M = 10 15 g/s. In other words, if the accretion rate is 
less than this critical value, a NS will accrete the mass 
less than 0.1 M Q and may evolve to a recycled PSR with 
a slightly longer spin period than 20 ms. Furthermore, 
we stress that this critical minimum accretion rate is an 
averaged value in the whole binary accretion history. 



(3) It is noticed that the position of a pulsar and 
its spin-up line in the B-P diagram are slightly affected 
by the choices of the NS mass and radius, and the ra- 
dius has more influence than the mass. However, un- 
like the measurements of the NS masses, the NS radius 
has no t yet been satisfactorily d etermined from observa- 
tions (IZhang et al.ll2007U201ll ). It is usually estimated 
to be in the range of 10-20 km, thus we cannot sure if 
the standard value of R = 10km introduces a big error 
in defining the PSR position in the B-P diagram. 

(4) The classification of the BPSRs by the conditions 
of the companion mass, eccentricity and orbit period 
have given us a clear illustrations of their properties in 
B-P diagram. Most millisecond BPSRs with small com- 
panion masses (< 0.45M Q ) should be long lived, such 
as 10 8_9 yr, while ~ 0.1 — 0.2M Q is accreted during the 
accretion phase. In addition, the majority millisecond 
BPSRs share the small eccentricities e < 0.01, which 
implies that their orbits were circularized during the 
recycle process. In fact, there have already been many 
researches on this, e.g. by Tauris & Sacinije (1999), 
Tauris, van den Heuvel & Savonijc (2000), and Tau- 
ris, Langer & Kromer (2012). However, there is no 
clear and simple relation between the orbital period 
and BPSR evolution, which hints that the evolution 
of the orbital period may be complicated in the recycle 
process. 

(5) To interpret why the seven BPSRs lie above the 
limit spin-up line, we start with the binary histories of 
these sources. The first four BPSRs have the heavy star 
companions (M c > 4M Q ) staying at the main sequence 
phases, and they are still spinning down via magnetic 
dipolc radiation. The recycle processes of these systems 
have not yet started, so they still go on the evolution- 
ary tracks of the first pulsar spin-down. The fifth BPSR 
lies in a double NS system (NS+NS). However, unlike 

— the system PSIL 1913+16 in which the recycled pul- 
— sai nasijeeii uo"rved, the observed pulsar J1906+0746 
should be the normal and the secondly formed pulsar. 
Namely, its companion should be a recycled NS that 
has not yet been observed as a pulsar. Thus this PSR 
is experiencing the spin-down by the electromagnetic 
radiation towards the death line. Both the sixth and 
seventh pulsars are similar, which have the degenerate 
white dwarfs (NS+WD) as companions. They might be 
the newly formed recycled pulsars with the Eddington 
rate, while either the particular accretion torque and 
disk structure or large NS radii may be responsible for 
their slightly passing through the spin-up line. 

In summary, the study of BPSRs in B-P diagram has 
provided us the clues on their accreting evolution his- 
tories, which is a helpful tool in constructing a link be- 
tween the radio MSPs and their progenitors, where the 



Table 1 Parameters of Seven Binary Pulsars above the Spin-up 


Line 






No. Name 


P/ms 


Porb/d 


e 


Mg/Mq 


TJyr 


B/G 


type 


1 B1259-63 


47.763 


1236.724 


0.8699 


4.14 


332000 


3.34E11 


NSMS 


2 J1638-4725 


763.933 


1940.9 


0.955 


8.08 


2.53E6 


1.93E12 


NSMS 


3 J0045-7319 


926.276 


51.1695 


0.8079 


5.27 


3.29E6 


2.06E12 


NSMS 


4 J1740-3052 


570.31 


231.0297 


0.5789 


15.82 


354000 


3.86E12 


NSMS 


5 J1906+0746 


144.072 


0.166 


0.0853 


1.37 


113000 


1.73E12 


DNS 


6 B1820-11 


279.829 


357.762 


0.7946 


0.78 


3.22E6 


6.29E11 


NSWD 


7 J1141-6545 


393.899 


0.1977 


0.1719 


1.02 


1.45E6 


1.32E12 


NSWD 



influences by their binary parameters can be properly 
investigated. The 186 binary pulsars from the 2008 pul- 
sar sample have been analyzed, and more samples are 
required to present robust conclusions on the evolutions 
of binary pulsars. With the recent construction of the 
powerful radio telescope, five hundred meter aperture 
spherical telescope (FAST), more than 4,000 new pul- 
sars a re expected to be observed in this decade (|Nan 
2006t ISmits et all I2009J: iNan et al.|[201ll) . while there 
will be a big sample set of binary pulsars, e.g. about 
500 ones, to show a clear trend of MSP evolution. 
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